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We study inverse statistical mechanics: how can one design a potential function so as to produce 
a specified ground state? In this paper, we show that unexpectedly simple potential functions suffice 
for certain symmetrical configurations, and we apply techniques from coding and information theory 
to provide mathematical proof that the ground state has been achieved. These potential functions 
are required to be decreasing and convex, which rules out the use of potential wells. Furthermore, 
we give an algorithm for constructing a potential function with a desired ground state. 
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I. INTRODUCTION 

How can one engineer conditions under which a desired 
structure will spontaneously self-assemble from simpler 
components? This inverse problem arises naturally in 
many fields, such as chemistry, materials science, biotech- 
nology, or nanotechnology (see for example ref. [IJ and 
the references cited therein) . A full solution remains dis- 
tant, but in this paper we develop connections with cod- 
ing and information theory, and we apply these connec- 
tions to give a detailed mathematical analysis of several 
fundamental cases. 

Our work is inspired by a series of papers by Rechts- 
man, Stillinger, and Torquato, in which they design 
potential functions that can produce a honeycomb Q, 
square cubic Q, or diamond [|| lattice. In this pa- 
per, we analyze finite analogues of these structures, and 
we show similar results for much simpler classes of po- 
tential functions. 

For an initial example, suppose 20 identical point par- 
ticles are confined to the surface of a unit sphere (in 
the spirit of the Thomson problem of how classical elec- 
trons arrange themselves on a spherical shell). We wish 
them to form a regular dodecahedron with 12 pentagonal 
facets. 

Suppose the only flexibility we have in designing the 
system is that we can specify an isotropic pair potential V 
between the points. In other words, the potential energy 
E\r{C) of a configuration (i.e., set of points) C is 
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In static equilibrium, the point configuration will assume 
a form that at least locally minimizes Ey(C). Can we 
arrange for the energy-minimizing configuration to be a 
dodecahedron? Furthermore, can we arrange for it to 
have a large basin of attraction under natural processes 
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FIG. 1: Two potential functions under which the regular do- 
decahedron minimizes energy: the blue one uses potential 
wells, while the green one is convex and decreasing. 



such as gradient descent? If so, then we can truly say 
that the dodecahedron automatically self-assembles out 
of randomly arranged points when the proper potential 
function is imposed. 

If we could choose V arbitrarily, then it would certainly 
be possible to make the dodecahedron the global mini- 
mum for energy by using potential wells, as in the blue 
graph in Fig. [T] By contrast, this cannot be done with 
familiar potential functions, such as inverse power laws, 
because the dodecahedron's pentagonal facets are highly 
unstable and prone to collapse into a triangulation. 

Unfortunately, the potential function shown in the blue 
graph in Fig.[T]is quite elaborate. Actually implementing 
precisely specified potential wells in a physical system 
would be an enormous challenge. Instead, one might ask 
for a simpler potential function, for example one that 
is decreasing and convex (corresponding to a repulsive, 



decaying force). 
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FIG. 2: The paths of points converging to the regular dodec- 
ahedron under the green potential function from Fig. [T] Only 
the front half of the sphere is shown. 

In fact, V can be chosen to be both decreasing and 
convex. The green graph in Fig.[T]shows such a potential 
function, which is described and analyzed in Theorem 2] 
We prove that the regular dodecahedron is the unique 
ground state for this system. We have been unable to 
prove anything about the basin of attraction, but com- 
puter simulations indicate that it is large (we performed 
1200 independent trials using random starting configu- 
rations, and all but 6 converged to the dodecahedron). 
For example, Fig. [2] shows the paths of the particles in 
a typical case, with the passage of time indicated by the 
transition from yellow to red. 

Our approach to this problem makes extensive use of 
linear programming. This enables us to give a probabilis- 
tic algorithm for inverse statistical mechanics. Using it, 
we construct simple potential functions with counterin- 
tuitive ground states. These states are analogues of those 
studied in refs. @, 0, 0, ll| , but we use much simpler po- 
tential functions. Finally, we make use of the linear pro- 
gramming bounds from coding theory to give rigorous 
mathematical proofs for some of our assertions. These 
bounds allow us to prove that the desired configurations 
are the true ground states of our potential functions. By 
contrast, previous results in this area were purely exper- 
imental and could not be rigorously analyzed. 



A. Assumptions and Model 

To arrive at a tractable problem, we make four funda- 
mental assumptions. First, we will deal with only finitely 
many particles confined to a bounded region of space. 
This is not an important restriction in itself, because pe- 
riodic boundary conditions could create an effectively in- 



finite number of particles. 

Second, we will use classical physics, rather than quan- 
tum mechanics. Our ideas are not intrinsically classical, 
but computational necessity forces our hand. Quantum 
systems are difficult to simulate classically (otherwise the 
field of quantum computing would not exist), and there 
is little point in attempting to design systems compu- 
tationally when we cannot even simulate them. Fortu- 
nately, classical approximations are often of real-world 
as well as theoretical value. For example, they are ex- 
cellent models for soft matter systems such as polymers 
and colloids [f| 0] ■ 

Third, we restrict our attention to a limited class of po- 
tential functions, namely isotropic pair potentials. These 
potentials depend only on the pairwise distances between 
the particles, with no directionality and no three-particle 
interactions; they are the simplest potential functions 
worthy of analysis. For example, the classical electric 
potential is of this sort. We expect that our methods 
will prove useful in more complex cases, but isotropic 
pair potentials have received the most attention in the 
literature and already present many challenges. 

Finally, we assume all the particles are identical. This 
assumption plays no algorithmic role and is made purely 
for the sake of convenience. The prettiest structures are 
often the most symmetrical, and using identical particles 
facilitates such symmetry. 

We must still specify the ambient space for the par- 
ticles. Three choices are particularly natural: we could 
study finite clusters of particles in Euclidean space, con- 
figurations in a flat torus (i.e., a region in space with 
periodic boundary conditions, so the number of particles 
is effectively infinite), or configurations on a sphere. Our 
algorithms apply to all three cases, but in this paper we 
will focus on spherical configurations. They are in many 
ways the most symmetrical case, and they are commonly 
analyzed, for example in Thomson's problem of arrang- 
ing classical electrons on a sphere. 

Thus, we will use the following model. Suppose we 
have N identical point particles confined to the surface 
of the unit sphere S" i_1 = {x € 1" : \x\ = 1} in n- 
dimensional Euclidean space R n . (We choose to work in 
units in which the radius is 1, but of course any other 
radius can be achieved by a simple rcscaling.) We use a 
potential function V : (0, 2] — * R, for which we define the 
energy of a configuration C C S n ' 1 as in Eq. (JTJ) . We 
define V only on (0, 2] because no other distances occur 
between distinct points on the unit sphere. 

Note that we have formulated the problem in an ar- 
bitrary number n of dimensions. It might seem that 
n = 3 would be the most relevant for the real world, 
but n = 4 is also a contender, because the surface of the 
unit sphere in R is itself a three-dimensional manifold 
(merely embedded in four dimensions). We can think of 
S 3 as an idealized model of a curved three-dimensional 
space. This curvature is important for the problem of 
"geometrical frustration" Q: many beautiful local con- 
figurations of particles do not extend to global config- 
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urations, but once the ambient space is given a small 
amount of curvature they piece together cleanly. As the 
curvature tends to zero (equivalently, as the radius of the 
sphere or the number of particles tends to infinity), we re- 
cover the Euclidean behavior. Although this may sound 
like an abstract trick, it sometimes provides a strikingly 
appropriate model for a real-world phenomenon; see, for 
example, figure 2.6 in ref. [§J, which compares the ra- 
dial distribution function obtained by X-ray diffraction 
on amorphous iron to that from a regular polytope in S 3 
and finds an excellent match between the peaks. 

The case of the ordinary sphere S 2 in R 3 is also more 
closely connected to actual applications that it might 
at first appear. One scenario is a Pickering emulsion, 
in which colloidal particles adsorb onto small droplets 
in the emulsion. The particles are essentially confined 
to the surface of the sphere and can interact with each 
other, for example, via a screened Coulomb potential or 
by more elaborate potentials. This approach has in fact 
been used in practice to fabricate colloidosomes @. See 
also the review article [l(|, in particular section 1.2, and 
the references cited therein for more examples of physics 
on curved, two-dimensional surfaces, such as amphiphilic 
membranes or viral capsids. 

B. Questions and Problems 

From the static perspective, we wish to understand 
what the ground state is (i.e., which configuration mini- 
mizes energy) and what other local minima exist. From 
the dynamic perspective, we wish to understand the 
movement of particles and the basins of attraction of the 
local minima for energy. There are several fundamental 
questions: 

1 . Given a configuration C , can one choose an isotropic 
pair potential V under which C is the unique ground 
state for \C\ points? 

2. How simple can V be? Can it be decreasing? Con- 
vex? 

3. How large can the basin of attraction be made? 

In this paper, we give a complete answer to the first 
question, giving necessary and sufficient conditions for 
such a potential to exist. The second question is more 
subtle, but we rigorously answer it for several important 
cases. In particular, we show that one can often use re- 
markably simple potential functions. The third question 
is the most subtle of all, and there is little hope of provid- 
ing rigorous proofs; instead, experimental evidence must 
suffice. 

The second question is particularly relevant for exper- 
imental work, for example with colloids, because only a 
limited range of potentials can be manipulated in the lab. 
Inverse statistical mechanics with simple potential func- 
tions was therefore raised as a challenge for future work 
in ref. [l[ . 




FIG. 3: The hypercube, drawn in four-point perspective. 

We will focus on four especially noteworthy structures: 

1. The 8 vertices of a cube, with 6 square facets. 

2. The 20 vertices of a regular dodecahedron, with 12 
pentagonal facets. 

3. The 16 vertices of a hypercube in four dimensions, 
with 8 cubic facets (see Fig. [3]). 

4. The 600 vertices of a regular 120-cell in four dimen- 
sions, with 120 dodecahedral facets (see Fig. 2]). 

The latter two configurations exist in four dimensions, 
but as discussed above the sphere containing them is a 
three-dimensional space, so they are intrinsically three- 
dimensional. 

These configurations are important test cases, because 
they are elegant and symmetrical yet at the same time 
not at all easy to build. The problem is that their facets 
are too large, which makes them highly unstable. Un- 
der ordinary potential functions, such as inverse power 
laws, these configurations are never even local minima, 
let alone global minima. In the case of the cube, one can 
typically improve it by rotating two opposite facets so 
they are no longer aligned. That lowers the energy, and 
indeed the global minimum appears to be the antiprism 
arrived at via a 45° rotation (and subsequent adjustment 
of the edge lengths). It might appear that this process 
always works, and that the cube can never minimize a 
convex, decreasing potential function. However, careful 
calculation shows that this argument is mistaken, and 
we will exhibit an explicit convex, decreasing potential 
function for which the cube is provably the unique global 
minimum. 
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FIG. 4: The Schlegel diagram for the regular 120-cell (with a 
dodecahedral facet in red). 

One reason why the four configurations mentioned 
above are interesting is that they are spherical analogues 
of the honeycomb and diamond packings from M 2 and M 3 , 
respectively. In each of our four cases, the nearest neigh- 
bors of any point form the vertices of a regular spherical 
simplex. They have the smallest possible coordination 
numbers that can occur in locally jammed packings (see, 
for example, ref. (Tlj|). 

II. POTENTIAL WELLS 

Suppose C is a configuration in S" 1 " 1 . The obvious 
way to build C is to use deep potential wells (i.e., lo- 
cal minima in V) corresponding to the distances between 
points in C, so that configurations that use only those dis- 
tances are energetically favored. This method produces 
complicated potential functions, which may be difficult 
to produce in the real world, but it is systematic and 
straightforward. In this section we rigorously analyze the 
limitations of this method and determine exactly when it 
works, thereby answering a question raised towards the 
end of ref. @ • 

The first limitation is obvious. Define the distance dis- 
tribution of C to be the function d such that d(r) is the 
number of pairs of points in C at distance r. The distance 
distribution determines the potential energy via 

*V(C) = £n?-Mr). (2) 

r 

Thus, C cannot possibly be the unique ground state un- 
less it is the only configuration with its distance distri- 
bution. 



The second limitation is more subtle. The formula ([2]) 
shows that Ey depends linearly on d. If d is a weighted 
average of the distance distributions of some other config- 
urations, then the energy of C will be the same weighted 
average of the other configurations' energies. In that 
case, one of those configurations must have energy at 
least as low as that of C. Call d extremal if it is an ex- 
treme point of the convex hull of the space of all distance 
distributions of |C|-point configurations in S 1 " -1 (i.e., it 
cannot be written as a weighted average of other distance 
distributions). If C is the unique ground state for some 
isotropic pair potential, then d must be extremal. 

For an example, consider three-point configurations on 
the circle S 1 , specified by the angles between the points 
(the shorter way around the circle). The distance dis- 
tribution of the configuration with angles 45°, 90°, and 
135° is not extremal, because it is the average of those 
for the 45°, 45°, 90° and 90°, f35°, 135° configurations. 

Theorem 1. If C is the unique configuration in S n 
with its distance distribution d and d is extremal, then 
there exists a smooth potential function V : (0,2] — > R 
under which C is the unique ground state among config- 
urations of \C\ points in S n . 

The analogue of Theorem Q] for finite clusters of parti- 
cles in Euclidean space is also true, with almost exactly 
the same proof. 

Proof: Because d is extremal, there exists a function I 
defined on the support supp(d) of d [i.e., the set of all r 
such that d(r) ^ 0] such that d is the unique minimum 
of 

t^^€(r)*(r) 

r 

among |C|-point distance distributions t with supp(f) C 
supp(d). Such a function corresponds to a supporting hy- 
perplane for the convex hull of the distance distributions 
with support contained in supp(d). 

For each e > 0, choose any smooth potential function 
V e such that V £ (r) — £(r) for r £ supp(d) while 

V e (s)>J2t(r)d(r) 

r 

whenever s is not within s of a point in supp(d). This 
is easily achieved using deep potential wells, and it guar- 
antees that no configuration can minimize energy unless 
every distance occurring in it is within £ of a distance 
occurring in C. Furthermore, when e is sufficiently small 
[specifically, less than half the distance between the clos- 
est two points in supp(d)], choose V e so that for each 
r S supp(d), we have V(s) > V(r) whenever \s — r\ < e 
and s^r. 

For a given e, there is no immediate guarantee that C 
will be the ground state. However, consider what hap- 
pens to the ground states under V e as e tends to 0. All 
subsequential limits of their distance distributions must 
be distance distributions with support in supp(<i). Be- 
cause of the choice of £, the only possibility is that they 



5 



are all d. In other words, as e tends to the distance 
distributions of all ground states must approach d. Be- 
cause the number of points at each given distance is an 
integer, it follows that when e is sufficiently small, for 
each r € supp(d), there are exactly d(r) distances in 
each ground state that are within e of r. Because V £ 
has a strict local minimum at each point in supp(d), it 
follows that it is minimized at d (and only at d) when 
s is sufficiently small. The conclusion of the theorem 
then follows from our assumption that C is the unique 
configuration in S" -1 with distance distribution d. □ 
It would be interesting to have a version of Theorem [1] 
for infinite collections of particles in Euclidean space, but 
there are technical obstacles. Having infinitely many dis- 
tances between particles makes the analysis more com- 
plicated, and one particular difficulty is what happens if 
the set of distances has an accumulation point or is even 
dense (for example in the case of a disordered packing). 
In such a case there seems to be no simple way to use 
potential wells, but in fact a continuous function with a 
fractal structure can have a dense set of strict local min- 
ima, and perhaps it could in theory serve as a potential 
function. 



III. SIMULATION-GUIDED OPTIMIZATION 

In this section we describe an algorithm for optimiz- 
ing the potential function to create a specified ground 
state. Our algorithm is similar to, and inspired by, the 
zero-temperature optimization procedure introduced in 
ref. Q ; the key difference is that their algorithm is based 
on a fixed list of competing configurations and uses simu- 
lated annealing, whereas ours dynamically updates that 
list and uses linear programming. (We also omit certain 
conditions on the phonon spectrum that ensure mechani- 
cal stability. In our algorithm, they appear to be implied 
automatically once the list of competitors is sufficiently 
large.) 

Suppose the allowed potential functions are the linear 
combinations of a finite set V\ , . . . , Vk of specified func- 
tions. In practice, this may model a situation in which 
only certain potential functions are physically realizable, 
with relative strengths that can be adjusted within a 
specified range, while in theory we may choose the basic 
potential functions so that their linear combinations can 
approximate any reasonable function arbitrarily closely 
as k becomes large. 

Given a configuration C C S*™ -1 , we wish to choose a 
linear combination 

v = \ 1 v 1 + --- + x k v k 

so that C is the global minimum for Ey. We may also 
wish to impose other conditions on V, such as monotonic- 
ity or convexity. We assume that all additional conditions 
are given by finitely many linear inequalities in the coef- 
ficients Ai, . . . , Afc. (For conditions such as monotonicity 
or convexity, which apply over the entire interval (0, 2] of 



distances, we approximate them by imposing these con- 
ditions on a large but finite subset of the interval.) 

Given a finite set of competitors Ci, .. . , Ct to C, we 
can choose the coefficients by solving a linear program. 
Specifically, we add an additional variable A and impose 
the constraints 

E v (Ci) >E V {C)+A 

for 1 < i < t, in addition to any additional constraints (as 
in the previous paragraph). We then choose Ax, • ■ • , Afe 
and A so as to maximize A subject to these constraints. 
Because this maximization problem is a linear program, 
its solution is easily found. 

If the coefficients can be chosen so that C is the global 
minimum, then A will be positive and this procedure 
will produce a potential function for which C has en- 
ergy less than each of C\, . . . ,Ce- The difficulty is how 
to choose these competitors. In some cases, it is easy to 
guess the best choices: for example, the natural competi- 
tors to a cube are the square antiprisms. In others, it 
is far from easy. Which configurations compete with the 
regular 120-cell in S 3 ? 

Our simulation-guided algorithm iteratively builds a 
list of competitors and an improved potential function. 
We start with any choice of coefficients, say Ai = 1 and 
Aa = • • • = Afc = 0, and the empty list of competitors. 
We then choose \C\ random points on S"™ -1 and minimize 
energy by gradient descent to produce a competitor to C, 
which we add to the list (if it is different from C) and 
use to update the choice of coefficients. This alternation 
between gradient descent and linear programming con- 
tinues until either we are satisfied that C is the global 
minimum of the potential function, or we find a list of 
competitors for which linear programming shows that A 
must be negative (in which case no choice of coefficients 
makes C the ground state) . 

This procedure is only a heuristic algorithm. When 
A is negative, it proves that C cannot be the ground 
state (using linear combinations of V\ , . . . , V k satisfying 
the desired constraints), but otherwise nothing is proved. 
As the number of iterations grows large, the algorithm 
is almost certain to make C the ground state if that is 
possible, because eventually all possible competitors will 
be located. However, we have no bounds on the rate at 
which this occurs. 

We hope that C will not only be the ground state, but 
will also have a large basin of attraction under gradi- 
ent descent. Maximizing the energy difference A seems 
to be a reasonable approach, but other criteria may do 
even better. In practice, simulation-guided optimization 
does not always produce a large basin of attraction, even 
when one is theoretically possible. Sometimes it helps to 
remove the first handful of competitors from the list once 
the algorithm has progressed far enough. 
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IV. RIGOROUS ANALYSIS 

The numerical method described in the previous sec- 
tion appears to work well, but it is not supported by 
rigorous proofs. In this section we provide such proofs in 
several important cases. The key observation is that the 
conditions for proving a sharp bound in Proposition[3] be- 
low are themselves linear and can be added as constraints 
in the simulation-guided optimization. While this does 
not always lead to a solution, when it does the solution 
is provably optimal (and in fact no simulations are then 
needed). The prototypical example is the following the- 
orem: 

Theorem 2. Let the potential function V: (0,2] — * K be 
defined by 

1 1.13 0.523 
V ^) = ^- — + — ■ 

Then the cube is the unique global energy minimum 
among 8-point configurations on S 2 . The function V is 
decreasing and strictly convex. 

Theorem [2] is stated in terms of a specific potential 
function, but of course many others could be found using 
our algorithm. Furthermore, as discussed in the con- 
clusions below, the proof techniques are robust and any 
potential function sufficiently close to this one works. 

The potential function used in Theorem [5] is modeled 
after the Lennard- Jones potential. The simplest general- 
ization (namely a linear combination of two inverse power 
laws) cannot work here, but three inverse power laws suf- 
fice. The potential function in Theorem [5] is in fact de- 
creasing and convex on the entire right half-line, although 
only the values on (0, 2] are relevant to the problem at 
hand and the potential function could be extended in an 
arbitrary manner beyond that interval. 

To prove Theorem[2J we will apply linear programming 
bounds, in particular, Yudin's version for potential energy 
[12] ] . Let Pi denote the ith degree Gegenbauer polynomial 
for S n ~ 1 [i.e., with parameter (n — 3)/2, which we sup- 
press in our notation for simplicity] , normalized to have 
Pi(l) — 1. These are a family of orthogonal polynomials 
that arise naturally in the study of harmonic analysis on 
S n ~ 1 . The fundamental property they have is that for 
every finite configuration C C S n , 



E Pii(x,v))>0- 



x,y£C 

(Here, (x, y) denotes the inner product, or dot product, 
between x and y.) See section 2.2 of ref. [l3| for further 
background. 

The linear programming bound makes use of an aux- 
iliary function h to produce a lower bound on poten- 
tial energy. The function h will be a polynomial h(t) = 
^j i=0 oiiPi{t) with coefficients ao,...,ay > 0. It will 
also be required to satisfy h{t) < V{y/2 — 2t) for all 



t e [—1,1). Note that ^/2~~2~i is the Euclidean dis- 
tance between two unit vectors with inner product t, be- 
cause | a; — y\ 2 — \x\ 2 + \y\ 2 — 2(x, y) = 2 — 2(x, y) when 
\ x \ = \v\ = 1- We view h as a function of the inner prod- 
uct, and the previous inequality simply says that it is a 
lower bound for V. 

We say the configuration C is compatible with h if two 
conditions hold. The first is that h(t) = V(y/2 - 2t) 
whenever t is the inner product between two distinct 
points in C. The second is that whenever oti > with 
i > 0, we have ^2 X y( - c Pi{(x, y)) — 0. This equation holds 
if and only if for every z S S" 1 " 1 , YlxeC z )) = 0- 

(The subtle direction follows from Theorem 9.6.3 in 
ref. H3.) 

Proposition 3 (Yudin [12]). Given the hypotheses listed 
above for h, every N -point configuration in S n ~ l has V- 
potential energy at least (N 2 a^ — Nh(\))/2. If C is com- 
patible with h, then it is a global minimum for energy 
among all \C\-point configurations in and every 

such global minimum must be compatible with h. 

Proof: Let C C be any finite configuration with 

N points (not necessarily compatible with h). Then 

Ev{C) = \ £ V{^2-2{x,y)) 

>\ E MM) 



Nh(l) l ^ , N 



Nh(l) 1 



x,y£C 
d 



i=0 x,yeC 



> 



Nh(l) a 



? E MM) 



x,y£C 



N 2 a - Nh(l) 



The first inequality holds because h is a lower bound 
for V, and the second holds because all the Pi-sums are 
nonnegative (as are the coefficients a,). The lower bound 
for energy is attained by C if and only if both inequalities 
are tight, which holds if and only if C is compatible with 
h, as desired. □ 

Proof of Theorem^ It is straightforward to check that 
V is decreasing and strictly convex. To prove that the 
cube is the unique local minimum, we will use linear pro- 
gramming bounds. 

Let h be the unique polynomial of the form 

h(t) = a + axP x {t) + a 2 P 2 (t) + a 3 P 3 {t) + a 5 P 5 (t) 



P4 is missing) such that h(i) agrees with 
Tj to order 2 at t — ±1/3 and to order 1 at 



(note that 
V(V2 

t = — 1. These values of t are the inner products between 



7 



distinct points in the cube. One can easily compute the 
coefficients of h by solving linear equations and verify 
that they are all positive. Furthermore, it is straightfor- 
ward to check that h(t) < V{^2 - 2t) for all t € [-1, 1), 
with equality only for t 6 { — 1, —1/3, 1/3}. 

The cube is compatible with h, and to complete the 
proof all that remains is to show that it is the only 8- 
point configuration that is compatible with h. Every such 
configuration C can have only — 1, —1/3, and 1/3 as inner 
products between distinct points. For each y € C and 
1 < i < 3, 

xEC 

If there are N t points in C that have inner product t with 
y, then 

P<(-l)JV_i + Pi(-l/3)JV_ 1/3 + P(l/3)7V 1/3 + P(l) = 

for 1 < i < 3. These linear equations have the unique 
solution N-i = 1, N±i/3 = 3. 

In other words, not only is the complete distance distri- 
bution of C determined, but the distances from each point 
to the others are independent of which point is chosen. 
The remainder of the proof is straightforward. For each 
point in C, consider its three nearest neighbors. They 
must have inner product —1/3 with each other: no two 
can be antipodal to each other, and if any two were closer 
together than in a cube, then some other pair would be 
farther (which is impossible). Thus, the local configura- 
tion of neighbors is completely determined, and in this 
case, that determines the entire structure. □ 



Theorem 4. Let the potential function V : (0, 2] 
defined by 



VM _ n , rt 5, (t + l) 2 (t-l/3) 2 (t+l/3) 2 (t 2 -5/9) 2 
1 ' { + J + 6(1 - £) 2 

where t = 1 — r 2 /2. Then the regular dodecahedron is the 
unique global energy minimum among 20-point configu- 
rations on S 2 . The function V is decreasing and strictly 
convex. 

The proof is analogous to that of Theorem [21 except 
that we choose hit) = (1 + t) 5 . The proof of unique- 
ness works similarly. Note that the potential function 
used in Theorem 3] is physically unnatural. It does not 
seem worth carefully optimizing the form of this poten- 
tial function when it is already several steps away from 
real-world application. Instead, Theorems [4] through [6] 
should be viewed as plausibility arguments, which prove 
that there exists a convex, decreasing potential while al- 
lowing its form to be highly complicated. 

Theorem 5. Let the potential function V: (0, 2] — > K be 

defined by 

-13 + 7M + 5t 2 + 7t 3 + t 5 



V(r) 



120 

7(t+l)(t + l/2) 2 f 2 (t-l/2) 2 
120(1 -t) 



where t = 1 — r 2 /2. Then the hypercube is the unique 
global energy minimum among 16-point configurations on 
S 3 . The function V is decreasing and strictly convex. 

For the 120-cell, let q(t) be the monic polynomial 
whose roots are the inner product between distinct points 
in the 120-cell; in other words, 

q{t) = t(t + l)(4t + l)(4t - l)(4t + 3)(4t - 3) 
(2t + l)(2t- 1)(16< 2 - 5) 
(4< 2 + 2t-l)(4t 2 -2i-l) 
(16i 2 + 4<-l)(16i 2 -4t-l) 
(16i 2 + At - ll)(16i 2 -it- 11) 
(16i 2 + 12t + l)(16t 2 - 12t + 1) 
(16i 2 + 20i + 5)(16i 2 - 20i + 5)/2 50 . 

Let mi, ... , m 2 g be the integers 2, 4, 6, 8, 10, 14, 16, 
18, 22, 26, 28, 34, 38, 46, 1, 3, 5, 7, 9, 11, 13, 15, 17, 19, 
21, 23, 25, 27, 29 (in order), and let Ci,...,ci 7 be 1, 2/3, 
4/9, 1/4, 1/9, 1/20, 1/20, 1/15, 1/15, 9/200, 3/190, 0, 
7/900, 1/40, 1/35, 3/190, and 1/285. 



Theorem 6. Let the potential function V : (0, 2] 
defined by 



be 



V(r)=J2c i P i (t) + J2 



10 6 



10' 



At) 2 
i-t' 



where t — 1— r 2 /2. Then the regular 120-cell is the unique 
be global energy minimum among 600-point configurations 
on S 3 . The function V is decreasing and strictly convex. 



The proofs of Theorems [5] and [H] use the same tech- 
niques as before. The most elaborate case is the 120-cell, 
specifically the proof of uniqueness. The calculation of 
the coefficients N t , as in the proof of Theorem^ proceeds 
as before, except that the P%2 sum does not vanish (note 
that the coefficient en of Pi 2 in V is zero). Nevertheless, 
there are enough simultaneous equations to calculate the 
numbers N t . Straightforward case analysis suffices to 
show then that the four neighbors of each point form 
a regular tetrahedron, and the entire structure is deter- 
mined by that. 



V. CONCLUSIONS AND OPEN PROBLEMS 

In this article, we have shown that symmetrical config- 
urations can often be built by using surprisingly simple 
potential functions, and we have given a new algorithm 
to search for such potential functions. However, many 
open problems remain. 

One natural problem is to extend the linear program- 
ming bound analysis to Euclidean space. There is no 
conceptual barrier to this (section 9 of ref. [HI develops 
the necessary theory) , but there are technical difficulties 
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FIG. 5: Gaussian energy on the space of two-dimensional 
lattices (red means high energy) . 

that must be overcome if one is to give a rigorous proof 
that a ground state has been achieved. 

A second problem is to develop methods of analyzing 
the basin of attraction of a given configuration under 
gradient descent. We know of no rigorous bounds for the 
size of the basin. 

The review article [l[ raises the issue of robustness: 
Will a small perturbation in the potential function (due, 
for example, to experimental error) change the ground 
state? One can show that the potential in Theorem [2] 
is at least somewhat robust. Specifically, it follows from 
the same proof techniques that there exists an e > 
such that if the values of the potential function and its 



first two derivatives are changed by a factor of no more 
than 1 + e, then the ground state remains the same. It 
would be interesting to see how robust a potential func- 
tion one could construct in this case. The argument 
breaks down slightly for Theorems [4] through [6l but they 
can be slightly modified to make them robust. 

It would also be interesting to develop a clearer geo- 
metrical picture of energy minimization problems. For 
example, for Bravais lattices in the plane, the space of 
lattices can be naturally described by using hyperbolic 
geometry [see, for example, ref. [Hj], (pp 124-125)]. Fig. [5] 
shows a plot of potential energy for a Gaussian potential 
function, drawn by using the Poincare disk model of the 
hyperbolic plane. Each point corresponds to a lattice, 
and the color indicates energy (red is high). The local 
minima in yellow are copies of the triangular lattice; the 
different points correspond to different bases. The sad- 
dle points between them are square lattices, which can 
deform into triangular lattices in two different ways by 
shearing the square along either axis. The red points on 
the boundary show how the energy blows up as the lattice 
becomes degenerate. In more general energy minimiza- 
tion problems, we cannot expect to draw such pictures, 
but one could hope for a similarly complete analysis, with 
an exhaustive list of all critical points as well as a descrip- 
tion of how they are related to each other geometrically. 
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